NASA TECHNICAL MEMORANDUM 100494 
USAAVSCOM TECHNICAL MEMORANDUM 87-B-12 

STRAIN-ENERGY-RELEASE RATE ANALYSIS 

D~ 33 

OF THE END-NOTCHED FLEXURE SPECIMEN r 
USING THE FINITE-ELEMENT METHOD 


S. A. Salpekar, I. S. Raju, and T. K. O'Brien 


N88-1 37 35 

Unclas 
0111863 

NOVEMBER 1987 


r NASA -TH- 10 04 9 4) STRAIN- ENERGY-RELEASE BATE 

ANALYSIS OF THE* END-NOTCHED FLEXURE SPECIMEN 

I “LnSiH-w-m ™ JNASA) 

33 P Avail: NT1S HC &03/HF 


A01 


CSCE 20K 




lunsn 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665 



STRAIN-ENERGY-RELEASE RATE ANALYSIS OF THE END-NOTCHED 


FLEXURE SPECIMEN USING THE FINITE -ELEMENT METHOD 

S. A. Salpekar* , I. S. Raju*, and T. K. O'Brien** 

NASA Langley Research Center 
Hampton, Virginia 23665-5225 

SUMMARY 

Two-dimensional finite-element analysis of the end-notched flexure 
specimen was performed using 8 -node isoparametric, parabolic elements to 
evaluate compliance and mode II strain-energy-release rates, Gjj. The Gjj 
values were computed using two different techniques: the virtual crack- 

closure technique (VCCT) and the rate of change of compliance with crack 
length (compliance derivative method) . The analysis was performed for 
various crack- length- to -semi -span (a/L) ratios ranging from 0.2 to 0.9. 

Three material systems representing a wide range of material properties were 
analyzed. The compliance and strain-energy-release rates of the specimen 
calculated with the present finite -element analysis agree very well with 
beam theory equations including transverse shear. The Gjj values 
calculated using the compliance derivative method compared extremely well 
with those calculated using the virtual crack-closure technique. The Gj j 
values obtained by the compliance derivative method using the top or bottom 
beam deflections agreed closely with each other. The strain-energy-release 
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rates from a plane-stress analysis were higher than the plane-strain values 
by only a small percentage, indicating that either assumption may be used in 
the analysis. The Gj j values for one material system calculated from the 
finite-element analysis agreed with one solution in the literature and 
disagreed with the other solution in the literature. 

v INTRODUCTION 

With the growing use of laminated composites in aircraft structures, 
more attention is being paid to the failure modes of such composites. A 
major source of failure in such materials is delamination (also referred to 
as a crack in this paper) along ply interfaces. In order to better 
understand the behavior of laminated composites and to be able to 
characterize and predict failure, the interlaminar fracture toughness must 
be evaluated. 

ASTM has undertaken a round robin test program to evaluate various test 
methods for determining the mode II interlaminar fracture toughness of 
laminated composite materials. The end-notched flexure (ENF) test [1] is 
being considered for the determination of mode II interlaminar fracture 
toughness . 

The ENF test is a three-point bend test of a unidirectional laminate 
and is used to obtain the compliance and mode II strain-energy-release rate 
(Gji) of the test specimen. A theoretical estimation of the specimen 
compliance and Gji is useful in correlating and predicting the results of 
experiments involving different material properties and specimen 
configurations. For this reason, beam theory equations were proposed in the 
literature for the ENF specimen [1,2]. The accuracy of these equations 
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needs to be established by methods which do not use simplifying assumptions 
of beam theory. Attempts were made to analyze the ENF specimen by two- 
dimensional finite -element analysis [3-6]. However, the strain-energy- 
release rates in reference 3 disagreed with those in reference 4 for several 
values of (a/L). 

The purpose of this paper is, therefore, to analyze the ENF specimen by 
finite-element analysis using higher order elements, to determine the 
specimen compliance and the strain-energy-release rates, and to compare 
these results with the beam theory solutions. Two methods, virtual crack 
closure technique (VCCT) and compliance derivative technique, were used to 
calculate the Gjj values. The Gjj values calculated from these two 
methods are compared to each other and to those in references 3 and 4. The 
analyses were performed for three commonly used material systems for values 
of delamination length to semi-span (a/L) ratios ranging from 0.2 to 0.9. 

SYMBOLS 

crack length 
beam width 

beam compliance considering bending deformation 
beam compliance considering bending and shear 
deformation 

beam compliance from finite-element analysis E^]^,E 33 
Young's modulii in x- and y-directions , respectively 

mode II strain-energy-release rate from beam 
theory considering bending and shear deformation 
transverse shear modulus 
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mode II strain- energy- release rate from finite - 
element analysis 
beam half - thickness 
beam half- span 
line load 

deflection under a unit load calculated using 
finite-element analysis 

deflection under a unit load calculated using 
beam theory with shear deformation 
Poisson's ratio 

length of the side of a square element around the 
crack tip 

ANALYSIS 

Specimen Configuration, Loading, and Materials 
Figure 1 shows the configuration of the specimen. The specimen is a 
24-ply unidirectional laminate with a through-width delamination at mid- 
depth at one end of the specimen. As shown in the figure, a is the 
delamination length, L is the specimen half- span, b is the width, and 
2h is the thickness. The half-depth, h, was assumed to be 0.066925 in., 
and the beam half-span was taken equal to 1.5 in. A line load, P, is 
applied at the mid- span as shown in Figure 1. 

Three sets of material properties which are shown in Table 1 were 
considered in this study. Material 1 properties were used to compare the 
present results with the results of reference 3. Material 2 represents 
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typical properties for S2/SP250 glass/epoxy, and material 3 represents those 
of AS4/PEEK. 


Beam Theory Solutions 

Simple expressions for compliance and in the ENF test are given 

in the literature [1,2]. These expressions were obtained by using beam 
theory with and without shear deformation. 

The equation for beam compliance, considering only bending deformation, 
was given by Russell [1] . For a mid-plane delamination in a unidirectional 
composite beam, this equation reduces to [5] 

C BT - (l + 1.5(a/L) 3 y^4E n b(h/L) 3 ^ (1) 

This equation for C BT was modified by Carlsson et al. [2] to include the 
effect of transverse shear deformation 



( 2 ) 


The strain- energy- release rates are calculated by taking the derivative 
of the compliance expression 
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The Gji values corresponding to compliance Equations (1) and (2), 
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respectively, are 

- (9a 2 P 2 )/(l6E 11 b 2 h 3 ) (4) 

G n- G n( 1 + 0 - 2 ( E u''°i3)< h ^) 2 ) < 5 > 

The compliance and strain-energy-release rate values from these simple 
expressions will be compared with the finite -element results. 

Finite -Element Analysis 

Model . - A two-dimensional plane-strain analysis of the specimen in 
Figure 1 was performed because the width is large in comparison to the 
thickness. The difference between plane-strain and plane-stress assumptions 
will be discussed later in this paper. 

Figure 2 shows a finite -element model with 1699 nodes and 516 8 -node 
isoparametric, parabolic elements. A very fine mesh was used in the 
vicinity of the delamination tip, with square elements around the 
delamination tip. The same model was used to study the effect of (a/L) 
ratio on Gjj. This was accomplished by moving the support points and the 
central load point to achieve the desired (a/L) ratio. The (L/h) ratio was 
the same in all the analyses. A wide range of (a/L) ratios (0.2 to 0.9) and 
material properties (presented in Table 1) were considered. 

Due to transverse loading and an (L/h) ratio of 22.4, the ENF specimen 
behaves like a plate. Therefore, to examine if reduced integration [7] is 
needed, all the elements in the model were integrated with three -point- 
Gaussian quadrature in the two coordinate directions. Then, the analysis 
was repeated with two-point-Gaussian quadrature. The deflection under the 
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load and the Gjj values calculated by the two quadrature schemes differed 
by less than 0.1 percent. Therefore, in the rest of the analysis, the two- 
point-Gaussian quadrature was employed. 

Crack-Face Constraints . - Initially, the finite -element model was 
analyzed allowing nodes on either side of the delamination to deform without 
any restraint. The resulting deflections indicated that delamination faces 
had crossed into each other. Because this deformation is physically 
inadmissible, the nodes along the crack face must be restrained so that such 
crossover is prevented [3,4]. When such a constraint is prescribed, the 
delamination faces press against each other and frictional forces will 
develop along the delamination faces. Thus, in a rigorous analysis, 
friction needs to be considered. However, the coefficient of friction (/i) 
needs to be assumed because data on /i are not experimentally available. 
Finite -element analyses were performed in references 4 and 5 for various 
values of friction, /i, and showed that the effect of friction is to lower 
the Gjj values, in general, by a small amount. Therefore, in the present 
analysis, friction is neglected, i.e. , smooth delamination faces are 
assumed. 

In the present analysis, two methods that prevent crossover of the 
delamination face nodes were used. Figure 3 illustrates these two methods. 
In the first method, the nodes on either side along the entire delamination 
face were constrained to move the same vertical distance (i.e., to have the 
same v-displacement) , but were allowed to slide past each other freely 
(i.e., to have different u-displacements) . The nodes on the top and bottom 
crack face (e.g., g,g') are shown separated in the vertical direction in 
Figure 3 for clarity only. In the second method, the same constraints as in 
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the first method were used, but only on a small region of the crack face. 
Consider the node r on the crack face vertically above the right support 
(see Fig. 3). The nodes p, q, r, s, t belong to two elements, one on 
either side of node r. Only these nodes were subjected to multipoint 
constraints; the remaining nodes on the crack face were unconstrained. The 
reason for this prescription is that the load transfer from the upper half 
of the beam to the lower half takes place in a very small region around 
node r [5]. Thus, only a few nodes around node r need to be 
constrained. The differences in the solutions obtained using the two 
methods will be discussed later. 

Strain-Energv-Release Rates .- The strain-energy-release rates were 
calculated using Irwin's virtual crack-closure technique. The equations 
utilized were similar to those given in references 8 and 9. 


G_ -= - . /v rt - v..\ + F ./v - v ,\1 (6a) 

I 2A| yi( i £'J yj( m m'H 

G _ « - If .(\i 0 - u \ + F ./u - u (6b) 

II 2A xi y H 2 J xjl m m' j I 

where A is the element size, Fxi and F yi are the forces at node i in 
the x- and y-directions , and (u^ - u ^') and (v^ - v g') are the relative 
sliding and opening displacements at node jg , respectively (see Fig. 4). 
Forces at node j and relative displacements between m and m' are 
defined similarly. 

A compliance derivative method was also used to obtain the strain- 
energy- release rates. For each (a/L) ratio, the finite-element solution 
yields the deflection under the unit load (and, hence, the compliance) of 
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the specimen. Using these numerical compliance values for various (a/L) 
ratios, a curve fit was made to obtain the compliance of the specimen as a 
function of (a/L) as 


C - A x + A 2 (a/L) + A 3 (a/L ) 2 + . . . 


(7) 


where A^, A 2 , A 3 , etc., are numerical constants 
square curve fit. The strain- energy- release rate, 
differentiating Equation (7) and then substituting 


determined by a least- 
s'll » was obtained by 
in Equation (3) 


2 

_ P dC 
II “ 2bL d(a/L$ 
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The G 11 values by this procedure were obtained using three-, four-, and 
five-term approximations in Equation (7). Comparison of Gji values 
indicated a maximum difference of 1.34 percent between the three sets of 
results. Thus, the cubic polynomial was used in the curve fitting for 
compliance. 


RESULTS AND DISCUSSION 

In this section, various factors that affect the finite -element 
solution for the ENF specimen are presented first. Then, the compliance and 
the strain- energy- release rates obtained by the various methods presented 
earlier are discussed. Next, the present results are compared to the beam 
theory solutions and other finite-element solutions. 

Multipoint Constraints 
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The normalized displacements and Gu values (computed by VCCT) 
obtained for material 1 from the two methods of prescribing multipoint 
constraints were found to be identical for configurations ranging from very 
small to very large crack lengths. No overlapping of the crack faces 
occurred behind the crack tip when the second method was used, indicating 
the absence of mode I type deformation. Comparison of the two methods 
confirms that the load transfer from the upper half section of the beam to 
the lower half occurs in the immediate vicinity of node r in 
Figure 3. This was also confirmed by summing the vertical forces at 
locations p through t in the first method and comparing the sum to the 
support reaction. Similar conclusions were reached in reference 5. The 
multipoint constraints for the entire crack face (method 1) were used in the 
rest of this study. 


Convergence Study 

In order to examine the convergence of the finite-element solution, the 
analysis was first performed using square elements of size 0.0075 in. (A/h - 
0.112) around the delamination tip. The elements around the delamination 
tip were then subdivided to obtain element sizes with (A/h) values of 0.056 
and 0.028 in the two subsequent analyses. Note that this subdivision 
extends through the specimen thickness and length along narrow bands. 

Figure 5 shows the convergence of Gjj values with mesh refinement for an 
(a/L) ratio of 0.4. 

The Gjj values changed by 1.4 percent as (A/h) changed from 0.112 to 
0.056. For the next refinement, the Gjj values changed by 0.5 percent as 
(A/h) changed from 0.056 to 0.028. These results indicate that Gji values 
with A/h - 0.028 can be considered as accurate. Therefore, the model with 
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A/h - 0.028 was used in all subsequent analyses. Converged Gjj values at 
other a/L ratios are also presented in Table 2 for completeness. 

Specimen Compliance 

The deflection under a unit central concentrated load (point D in Fig. 
1) equals the compliance of the specimen. The compliance values predicted 
with the finite -element analyses and those from the shear deformation theory 
(Eq. (2)) are presented in Table 3. Figure 6 shows typical results for 
material 1. Excellent agreement is observed between the two sets of 
results . 


Strain-Energy-Release Rate 

In this section, a comparison of strain-energy-release rates, Gxx, 
from two different techniques (VCCT and compliance derivative technique) is 
presented. The finite -element solutions are also compared with the beam 
theory solutions and the results in the literature. Finally, the effect of 
plane-stress and plane-strain assumptions on Gjj is examined. 

The Gjj values calculated by the virtual crack-closure technique (Eq. 
(6)) and the compliance derivative technique (Eq. (8)) for various materials 
are given in Table 4, and typical results are shown in Figure 7 for material 
1. The Gjj values calculated by the two methods are in excellent 
agreement, showing the consistency of the finite-element results. 

Often in experimental investigations, instead of measuring the 
deflection at the point of load application, some investigators prefer to 
measure the deflection at a point on the load line but on the bottom surface 
of the specimen (point D' in Fig. 1). The deflections at point D' will be 
referred to herein as the bottom deflections. These deflections are also 
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included in Table 3. The compliance derivative method was also applied to 
the bottom deflections and the Gji results are included in Table 4. 
Excellent agreement between Gjj values from the compliance derivative 
technique using either the top or bottom deflections was obtained, 
indicating that either deflection can be used. 

Comparison with Beam Theory Solutions .- In order to assess the 
usefulness of the beam theory equations in predicting strain- energy release 
rate, finite -element results for Gjj were compared with the beam theory 
solutions , Equations (4) and (5). These values are included in Table 4. 
This table and Figure 8 show that the finite -element analysis agrees better 
with the beam theory with transverse shear than with the theory without 
shear deformation. The finite-element Gjj values, are higher than the 
shear deformation theory, but the maximum difference is about 7 percent for 
the range of (a/L) ratios and materials considered. Note that this maximum 
difference did not occur at the same (a/L) ratio for the various materials 
considered herein. Also, the beam theory neglecting shear deformation gave 
lower Gji values compared to the finite-element values with larger 
percentage differences occurring at low values of (a/L) . 

Comparison with Other Solutions from the Literature .- Several finite - 
element solutions [3-5] of the ENF specimen exist in the literature. 

The specimen configurations in references 3 and 4 were modeled using 4-node 
quadrilateral elements. In reference 4, the multipoint constraints were 
prescribed all along the crack line similar to those in the first method in 
this study. In reference 3, "nonlinear truss elements" were used to 
prescribe the constraints along the crack line. In order to have a common 
basis for comparison, the results shown in Figure 3 of reference 4 and 
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Figure 7 of reference 3 were normalized with the G ■ values from Equation 

(5) and are shown with the results of the present analysis in Figure 9. 

Note that in references 3 and 4 specimens with different parameters, 
(E11/G13) and L, were analyzed. To be able to compare with the results of 
these references , the present analyses were repeated with the appropriate 
values of (E11/G13) and semi-span L. A new finite -element model for L - 2 
in. was developed. This model is very similar to that shown in Figure 2 
for L - 1.5 in. 

The results of the present analysis are in very good agreement with the 
results of Mall and Kochhar [4] and are in considerable disagreement with 
the results of reference 3. The disagreement between the present results 
and reference 3 results is larger for (a/L) values greater than 0.4. Also, 
the present results for the two sets of parameters showed very small 
differences . 

Gillespie, et al. [5], presented finite-element results for various 
values of (Eix/G^) and a semi-span, L - 1.5 in. The Gji values 
presented in reference 5 differed from the shear deformation beam theory by 
about 2 to 48 percent for (E11/G13) ratios ranging from 12.8 to 26.9. In 
contrast, in the present analysis, the maximum difference between the 
finite -element and the shear deformation beam theory is only about 7 
percent for ratios ranging from 10.5 to 31.7. The reasons for 

these discrepancies are not known. 

Plane Stress Versus Plane Strain . - The finite-element analysis in 
reference 3 was performed assuming plane-stress condition, whereas plane- 
strain condition was assumed in reference 4. To examine the differences due 
to these two assumptions, consider the expression for strain- energy- release 
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rate in Table 6 of reference 10 


G 



( 9 ) 


where 


c 



2A 12 + A 66 


1/2 


2A 
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( 10 ) 


for an orthotropic material. The Ajj's in Equation (10) are the compliance 
coefficients of the material. The stress- intensity factor Kjj is a 
function of specimen configuration and loading and is identical in plane - 
stress and plane-strain conditions. 

The ratio of Gjj for the plane -stress to the plane -strain case is 
therefore given by the corresponding ratio of c. This ratio was calculated 
to be 1.023 for S2/SP250 glass/epoxy (material 2) and 1.016 for AS4/PEEK 
(material 3). The corresponding ratios from finite-element analysis were 
obtained as 1.025 and 1.01, respectively. 

Since the difference in the plane stress and plane strain results are 
so small, either assumption may be used in an analysis. 


CONCLUDING REMARKS 

Two-dimensional finite-element analysis of the end-notched flexure 
specimen was performed using 8 -node isoparametric, parabolic elements to 
evaluate compliance and mode II strain- energy- release rates, Gjx- The Gn 
values were computed using two different techniques: the virtual crack- 

closure technique (VCCT) and the compliance derivative technique. The 
analysis was performed for various crack-length-to-semi-span (a/L) ratios 
ranging from 0.2 to 0.9. Three material systems covering a wide range of 
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material properties were considered. The compliance and Gjj values 
calculated from the finite-element analysis were compared with those based 
on beam theory equations and with other finite -element solutions from the 
literature. Based on this analysis, the following conclusions were reached: 

1. The compliance of the specimen calculated with the present finite- 
element analysis agrees very well (within 2.7 percent) with beam 
theory equations which consider shear deformation. The beam 
theory with shear deformation underestimates the Gjj values in 
comparison to finite -element results. The maximum difference was 
about 7 percent. Thus, beam theory equations which include the 
effect of shear deformation may be used to estimate the 
compliance and strain-energy-release rates of the end-notched 
flexure specimen. The results of beam theory equations which do 
not include the effect of shear deformation differ considerably 
from the finite -element results, especially at low (a/L) ratios. 

2. The Gjj values were also obtained from the compliance derivative 
method. The Gji values agreed extremely well (within 1 percent) 
with those calculated using the virtual crack- closure technique. 

3. A close agreement (within 0.3 percent) was obtained between Gjx 
results from the top or the bottom deflections by the compliance 
derivative method, indicating that either deflection may be used 
to measure compliance. 

4. The strain-energy-release rates from plane-stress analysis are 
higher than the plane-strain values by about 2.5 percent for 
S2/SP250 glass/epoxy and 1 percent for AS4/PEEK 
graphite/thermoplastic . 
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Table 1 - Material Properties Used in the Present Analysis 



E 

11 

MSI 

E 

33 

MSI 

G 

13 

MSI 

V 

13 

Material 1 

16.7 

1.4 

0.65 

0.3 

Material 2 

6.31 

2.5 

0.6 

0.25 

Material 3 

21.23 

1.5 

0.67 

0.37 
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Table 2 - Convergence of Gji Values with Mesh 
Refinement: (Material 1) 


a/L 


G x 10 6 in- lb/in 

> 

A/h 

- 0.1121 

0.0560 

0.0280 

0.2 


12.9415 

13.1274 

13.1970 

0.4 


45.4101 

45.9955 

46.2033 

0.6 


97.6077 

98.8145 

99.2334 

0.8 


169.368 

171.416 

172.118 

0.9 


211.564 

214.094 

214.959 
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Table 3 - Comparison of Compliance Values from Finite -Element 
Analyses and Beam Theory Solutions 



FE 


Beam Theory 

Percent 


Analysis 

Solution 

Difference* 

a/L 

x 10 5 

in- lb 

x 105 in- lb 

C FE . C SH 

xlOO 


Top of 
Beam 

Bottom of 
Beam 

** *** 
BT SH 

C FE 


Material 1 


0.0 

17.8895 

17.7610 

16.8551 

17.8895 

0.0 

0.2 

18.2047 

18.0734 

17.0574 

18.2470 

-0.23 

0.4 

19.8769 

19.7484 

18.4732 

19.8180 

0.30 

0.6 

24 . 1446 

24.0173 

22.3162 

23.8162 

1.36 

0.8 

32.2045 

32.0777 

29.7999 

31.4550 

2.33 

0.9 

38.0182 

37.8909 

35.2862 

37.0189 

2.63 


Material 2 


0.0 

44.7062 

44.5981 

44.6087 

45.7293 

-2.29 

0.2 

45.3901 

45.2820 

45 . 1440 

46.4327 

-2.30 

0.4 

49.4907 

49.3866 

48.8911 

50.3480 

-1.73 

0.6 

60.1585 

60.0563 

59.0619 

60.6868 

-0.88 

0.8 

80.5224 

80.4211 

78.8681 

80.6612 

-0.17 

0.9 

95.3141 

95.2129 

93.3883 

95.2654 

0.05 
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Table 3 - ( Continued ) . 



FE 

Analysis 

Beam Theory 
Solution 

Percent 

Difference* 

a/L 

x 10 5 in- lb 

x 10 5 in- lb 

C FE . C SH 

xlOO 


Top of Bottom of 

Beam Beam 

-kick 

BT SH 

C FE 

Material 3 


0.0 

14.2584 

14.1369 

13.2586 

14.2622 

-0.03 

0.2 

14.5045 

14.3830 

13.4177 

14.5718 

-0.46 

0.4 

15 . 8404 

15.7217 

14.5315 

15.8361 

0.03 

0.6 

19.2291 

19.1116 

17.5544 

19.0096 

1.14 

0.8 

25.6082 

25.4913 

23.4413 

25.0470 

2.19 

0.9 

30.1973 

30.080 

27.7569 

29.4379 

2.52 


C is compliance determined from displacement at the top of the beam. 
BT is beam theory considering bending deformation. 

SH is beam theory considering bending and shear deformation. 
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Table 4 - Comparison of Gjj Values Calculated by Various Methods 


FE Analysis Beam Theory Solution Percent 

x 10^ in-lb/in^ x 10^ in-lb/in^ Difference* 


FE SH 
G - G 

Compliance Compliance ** *** II II 

a/L VCCT Derivative Derivative BT SH x 100 

FE 

Top Bottom G 


II 


Material 1 


0.2 

13.1970 

13.2381 

13.2399 

10.113 

12.699 

3.77 

0.4 

46.2033 

46.2042 

46.2407 

40.452 

43.038 

6.85 

0.6 

99.2334 

99.2827 

99.3130 

91.017 

93.6038 

5.67 

0.8 

172.118 

172.474 

172.457 

161.809 

164.395 

4.49 

0.9 

214.959 

216.611 

216.556 

204.789 

207.376 

3.53 


Material 2 


0.2 

31.2631 

31.1766 

31.2102 

26.7652 

29.5668 

5.43 

0.4 

114.422 

114.386 

114.435 

107.061 

109.862 

3.99 

0.6 

249.717 

249 . 846 

249.884 

240.887 

243.688 

2.41 

0.8 

437.108 

437.559 

437.555 

428.243 

431.045 

1.39 

0.9 

549.781 

551.010 

550.975 

541.995 

544.797 

0.91 
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Table 4 - ( Continued ) . 



FE Analysis 

Beam Theory Solution 

Percent 


x 10^ in- lb/in 2 

X 10 6 

in- lb/in 2 

Difference* 

a/L 

Compliance Compliance 
VCCT Derivative Derivative 

** 

BT 

*** 

SH 

FE SH 

G - G 

II II 

x 100 

Top Bottom 



FE 

G 


II 


Material 3 


0.2 

10.6385 

10.5896 

10.6133 

7.9552 

10.4641 

1.64 

0.4 

36.7811 

36.8115 

36.8468 

31.8207 

34.3296 

6.67 

0.6 

78.6305 

78.7069 

78.7311 

71.5966 

74.1055 

5.76 

0.8 

135.966 

136.276 

136.266 

127.283 

129.792 

4.54 

0.9 

169.350 

170.938 

170.904 

161.092 

163.601 

3.40 


* FE 
G II 


** 


*** 


BT 

SH 


used here is obtained using VCCT. 

is beam theory considering bending deformation. 

is beam theory considering bending and shear deformation. 


i 

i 
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Fig. 1: Specimen configuration and loading. 
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Fig. 2: Finite-element model. 


Constraints: 
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Fig. 3: Effect of multipoint constraints. 







28 


and shear deformation beam theory (material 1). 



29 


VCCT and compliance derivative method (material 1). 



X 10 





i 
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Fig. 8: Comparison of Gy calculated using finite elements 
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Fig. 9: Comparison of Gy from various analyses. 
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